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ABSTRACT 

We investigate the implications of our measurement of the Lyman-a forest opacity at redshifts 2 < 
z < 4.2 from a sample of 86 high-resolution quasar spectra for the evolution of the cosmic ultraviolet 
luminosity density and its sources. The derived hydrogen photoionization rate F is remarkably flat 
over this redshift range, implying an increasing comoving ionizing emissivity with redshift. Because 
the quasar luminosity function is strongly peaked near z ~ 2, star-forming galaxies likely dominate 
the ionizing emissivity at z > 3. Our measurement argues against a star formation rate density 
declining beyond z 3, in contrast with existing state-of-the-art determinations of the cosmic star 
formation history from direct galaxy counts. Stellar emission from galaxies therefore likely reionized 
the Universe. 

Subject headings: Cosmology: diffuse radiation — methods: data analysis — galaxies: formation, 
evolution, high-redshift — quasars: absorption lines 
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1. INTRODUCTION 

The opacity of the Lyman-a (Lya) forest is set by 
a competition b etween hydrogen photo ionizations and 
recombinations (jCunn fc Petersonl Il965l ) and can thus 
serve as a direct probe of the photonization rate (e.g. 
iRauch et al.lll997l ). The hydrogen photoionization rate 
F is a particularly valuable quantity as it is an integral 
over all sources of ultraviolet (UV) radiation in the Uni- 
verse, 



F(z) = 47r r ^ J. 



(z)ct(z/). 



(1) 



where Ji, is the angle-averaged specific intensity of the 
background, cr(j^) is the photoionization cross section of 
hydrogen, and the integral is from the Lyman limit to 
infinity. As such, it bears a signature of cosmic stellar 
and quasistellar activity that is not subject to the com- 
pleteness issues to which direct source counts are prone. 
Moreover, unlike the redshifted radiation backgrounds 
observed on Earth, the Lya forest is a local probe of the 
high-redshift UV radiation, as only sources at approxi- 
mately the same redshift co ntribute to F at any point in 
the forest (e.g., Haardt & Madau 1991). In addition to 
being a powerful probe of galaxy formation and evolu- 
tion and a fu ndamental ingred ient of cosmological sim- 
ulations (e.g.. lEfstathiouiri992f ). identifying the sources 
that contribute most to the UV background is key to our 
understanding of the rcionization history of the Universe. 

In this Letter, we derive the photoionization rate im- 
plied by our measurement of the Lya forest opacity at 
2 < z < 4.2 from a sample of 86 high -resolution quasar 
spectra (jFaucher-Giguere et al.l[2008a ). for the first time 
consistently analyzing such a large data set (corrected 
for both continuum bias and metal absorption) over this 
redshift interval. We discuss the implications of its flat- 
ness over this redshift range for the relative contribu- 
tion of quasars and star-forming galaxies to the high- 
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redshift cosmic UV back ground. Throughout, we assume 
a WMAP5 cosmology (jKomatsu et al.l 120081 ). The full 



details of our analysis, a s well as supporting argument s, 
are presented elsewhere (jFaucher-Giguere et al.ll2008al) . 

2. THE PHOTOIONIZATION RATE FROM THE 
Lya FOREST 

The specific measurement we use is that of the 
Lya effective optical depth Teg in Az = 0.2 bins cor- 
rected for continuum bias and for metal a bsorption fol- 
lowing the results of ISchave et al.l (j2003f ) reported by 
iFaucher-Giguere et al.l (j2008cD . 

The effective optical depth is defined as 

T,s^-\n[{F){z)l (2) 

where {F) = (exp (— r)) is the mean t ransmission of the 
forest at redshift z and r is the local iGunn fc Petersonl 
(jl965( ) optical depth. In photoionization equilibrium and 
for a power-law temperature-density relation for the low- 
density intergalactic medium (IGM) of the form T ~ 
To{l + 5f (jHui fc Gnedin|[T99l . 



(3) 



with 



A{z) 



^Lya 



X x(x + o.sr) 



-0.7 



F 



H{z) 
-{l + zf 



(4) 



Here, 
sition 



is the oscillator strength of the Lya tran- 

^yQ, is its frequency, X and Y are the mass 

fractions of hydrogen and heliurn (resp ectively taken 
to be 0.75 and 0.25; iBurles et all 120011 ). Rq = 4.2 x 
10-13 cm? s-V(104 K)-o ^ and To is the IGM temper- 
ature at mean density (5 = 0). This expression is valid 
when all the intergalactic helium is fully ionized; an error 
< 8% may arise prior to Hell reionization. 

Given a volume-weighted probability density function 
(PDF) for the gas density A = 1 + (5, 



{F){z) 



dAP(A;z)exp(-r) 



(5) 
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Fig. 1. — (A) Photoionizat ion rate F infer r ed fro m our Lya effective optical depth measurement (black squares). The contribution from 
quasars calculated using the IHopkins et al.l ||2007|) quasar luminosity function is shown by the short dashes. (B) Comoving UV specific 
emissivity at 1500 A obtained by integrating galaxy UV luminosity functions. The green points are from 'Sawicki & Thompson' f2003) 
(Keck Deep Fields), the cyan points fr om Rcddy ct al. (2008) (Keck LBG), the red points from Yoshida ct al. (2006) (Subaru Deep Field), 
the magenta points from lBouwens e~a l. (2007) (Hubble Ultra Deep Field and other deep HST fields), and the blue point is based on the 
ISteidel et alji (1999) z ~ 4 LBGs, with the characteristic magnitude and faint-end slope set to the z ~ 3 values of [Re dd v ct al. (200J5(). See 
the text for caveats about the error bars shown. The black points show the emissivity implied by our Lya forest measurement, where the 
normalization was set to match the LF-derived values. (C) Comoving star formation rat e densit y implied by the UV emissivity (see text 
for details). Same color scheme as in (B). The long dashed curve shows the best fit of iHopkins fc Beacomi (120061 ) to the star formation 
history. 



We use the an alytical fit to gas-dynamical simulations of 
iMiralda-EscudeiaLl HOOO) for this PDF. For the IGM 
temperature, we inter polate between the Tp ^ 2 x 10^ 
K values measured by IZaldarriaga et al.l (|2001h from the 
Lya forest power spectrum, and /3 = 0.62, appropriate in 
the l imit of early hydrogen reionization (iHui &: GnedinI 

mn). 

With the above, we solve for the unique F that re- 
produces the measured TeS at each redshift. The results 
are shown in Figure [TJ^, where F (0.5 ± 0.1) x 10~^^ 
s""'^ is seen to be remarkably flat over the redshift range 
2 < z < 4.2. Note, however, that the absolute nor- 
malization of F depends on the cosmology, the thermal 
history of the IGM, as well as on the gas density dis- 
tribution (e.g.. iBolton et all 120051 ). and that a signifi- 
cant scatter between the results of different studies em- 
ploying the same basic meth od remains (see Figure 1 of 
iFaucher-Giguere et al.l l2008a). These systematic sources 
of uncertainty are not included in our analysis. On the 
other hand, the measurement we present consistently 
samples a large redshift interval with independent statis- 
tical errors and should therefore reliably trace the red- 
shift evolution of F. 

3. THE IONIZING SOURCES 

The UV background is generally assumed to be pro- 
duced by quasars and star-forming galaxies, but the rel- 
ative importance of these two populations remains un- 
certain. Moreover, it is unclear whether all the sources 
responsible for the ionization state of the high-redshift 
IGM are presently accounted for by magnitude-limited 
surveys. 

In Figure [T]A., we show the contribution of quasars 
as calculated using the i3-band realization of the 
ins et al.l ()2007h bolometric quasar luminosity func- 
tion (LF). The curve, subject to overall normalization 
uncertainties in the mean free path of ionizing photons, 
the spectral energy distribution of quasars, and the frac- 



tion of ionizing photons that they emit that escape into 
the IGM, has been renormalized to approximately match 
the total photoionization rate of the Lya forest at z = 2. 
Its shape is however robustly constrained at redshifts 
z > 2, owing to both an increasing dominance of the 
brightest quasars to the UV background and obscura- 
tion corrections decreasing in importance with redshift. 
In particular, the quasar contribution to F is strongly 
peaked near z = 2 and even if these objects produce 
the entire ionizing background at this redshift, they fall 
short of accounting for the total F measured at z = 4 by 
a factor > 5. 

To estimate the contribution of star-forming galaxies 
to the UV background, we consider recent determina- 
tions of the galaxy UV LF from Lyman break galaxy 
(LBG) surve ys bv iSawicki fc ThompsonI (|2006l ) (Keck 
Deep Fields) . iBOTawens et al.l ()2007[ ) (Hubble Ultra Deep 
Field [HUDF] and oth er deep Hubble Sp ace Telescope 
[HST] deep fields), Steidel ejall (il999,) andfReddv et all 
mm (Keck LBG), and lYoshida et al.1 (|200i ) (Subaru 
Deep Field). These were selected to be the most up- 
to-date measurements in the fields covered. The LFs 
have an effective wavelength near 1500 A, and the spe- 
cific emissivity at this wavelength is simply obtained by 
extrapolating and integrating them down to zero lumi- 
nosity. The error bars we quote are propagated from 
those on the individual Schechter parameters; because 
the latter are generally correlated, these will overestimate 
the tru e errors on the luminosity d ensities. Exceptions 
are the ISawicki fc ThompsonI (|2006f l points, for which we 
take the total luminosity densities and errors reported by 
the authors. In order to compare the UV luminosity den- 
sity of LBGs with our measured F, we convert the latter 
to a comoving emissivity at 1500 A. 

Given the proper mean free path of ionizing photons 
at the Lyman limit Amfp, we can calculate the total co- 
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moving specific emissivity implied by the measured F: 



h{am + 3) 
£912 ~ — : — ^ U 



^-3T 



(6) 



O'HI'^mfp 

(e.g., ISchirber & Bullockl 120031 ) for an ionizing back- 
ground with a power-law spectrum Jhi oc Here, 
(Thi is the photoionization cross section of hydrogen at 
the Lyman limit. Letting auv be the spectral index be- 
tween 912 A and 1500 A, we can calculate the emissivity 
at the wavelength probed by the galaxy UV luminosity 
functions, 



£1500 



1500 A 
7^;^ V 912 A 



£912, 



(7) 



where the escape fraction /esc accounts for the discon- 
tinuity at the Lyman limit owing to Lyman-continuum 
absorption associated with the host galaxy. 

The exact value for /esc is not well constrained at 
prese nt, but is likely to be at m ost a fe w percent 
(e^ lSteidel et al.ll200ll : IShaplev et"a l. 2006; Chen et all 
I2OO7D . Here, we simply note that £1500 oc F(l -I- z) 
for a mean free path Amfp oc {1 + z)'^, as appropri- 
ate if the i ncidence of Lyman-lirnit sys tems increases as 
(1 + z)^-^ (|Stengler-Larrea et al.l 1199 5V This is a con- 
servative assumption, as our conclusions regarding the 
need for an increasing comoving ionizing emissivity with 
redshift would only be strengthened if the absorbers re- 
sponsible for limiting the mean free path instead evolve 
as (1 -I- z)^, as is often assumed on the basis of the better 
studi ed lower column density systems (e.g.. Mada u et all 
|1999|) . We then solve for the normalization that mini- 
mizes the difference between the emissivities calcu- 
lated from the UV luminosity functions and our Lya for- 
est me asurem ent. Th e result is shown Figure [TJ3. In 
iFauch er-Gigue re et al.l (|2008a[ ). we show that for fidu- 
cial assumptions, only a small /esc ^ 0.5% is required 
for LBGs to solely account for the z ~ 3 ionizing back- 
ground. 

Within the large scatter, the redshift evolution of UV 
emissivity derived from the Lya forest is reasonably re- 
produced by the emission from LBGs only. The only hint 
of a decline of the galaxy UV emiss ivity near z = A comes 
from t he highest-redshift point of ISawicki fc Thompson! 
()2006l ) . This measurement is inconsistent with the higher 
poin ts from the Su baru Deep Field ( Yosh ida et al.|[2006r ) 
and ISteidel et al.l (1999). This may owe to cosmic 
variance in the relatively small Keck Deep Fields (169 
arcmin^ vs. ~850 arcmin^ for Subaru and Steidel et al.), 
or perhaps to selection effects (I Stanwav et al.l [20081 ). 

Quasars being clearly insufficient to solely account for 
the entire ionizing background implies that galaxies al- 
most certainly dominate at z > 3. Measurements of 
HeH to HI column density ratios however suggest that 
quasars do contribute a significant, perhaps dominant, 
fracti on of the ionizing background at their z ^ 2 peak 
(e.g.. IFaucher-Giguere et al.ll2008aD . 

4. THE COSMIC STAR FORMATION HISTORY 

The above results have interesting implications for the 
cosmic star formation history. In particular, several au- 
thors have previously found evidence for a peak in the 
star formation rate (SFR) density n ear z ~ 2 — 3 (e.g., 
iMadau et all 119981 : iHopkinsI l2004t iHopkins fc BeacomI 



l2006f ). Barring redshift evolution of dust obscuration, 
the escape fraction, or the initial mass function (IMF) 
of stars, the SFR density should trace the UV emis- 
sivity, which the Lya forest suggests instead increases 
continuously from z = 2toz = 4.2. As a represen- 
tative example of state-of-the-art determinations of the 
hi gh-redshift star forrnation history, we consider the fit 
of IHopkins fc BeacomI (|2006D to a large compilation of 
galaxy surveys. 

In Figure [Tp, we show the comoving SFR density we 
derived from the UV emissivities from both the Lya for- 
est (assuming that it arises solely from galaxies) and from 
direct measurements of the galaxy UV LF. We convert 
from specific UV emissivity to SFR density using 



= 1.08 X 10-28e,5oo, 



(8) 



where p^, is in units of comoving Mq yr"-'^ Mpc""^ 
provided £1500 is expressed in comoving erg s~^ Hz~^ 
Mpc~'^. This conversion is appropri ate for a "modified 
Salpet er A" IMF, consistent with the Hopkins fc BeacomI 
(|2006!) fit also shown on the Figure. Other IMFs 
would result in different conversion factors. However 
all the data points (and fit) in this plot would be 
equally renormalized and conclusions with regards to dis- 
crepant redshift evolutions would be unaffected. We ap- 
ply a UV obscuration correction factor of 3.4 over the 
entire redshift range, correspondi ng to the "common" 
obscuration correct i on ap plied by iHopkinsI ()2004D and 
IHopkins fc BeacomI ()2006[ ) at z > 3. Although this cor- 
rection is unlikely to be exact, it allows for a consistent 
comparison with the high-redshift fits to the star forma- 
tion history by these authors. 

We find no compelling evidence for a decline in the 
comoving SFR density over the redshift range probed 
by our measurement, either from it or from the di- 
rectly measured UV LF, in contrast to the best fit of 
ins fc BeacomI (|2006l ). Inspection of Figure [T]C sug- 
gests that the present data are instead roughly con- 
sistent with a constant ^ 0.2 Mq yr~^ Mpc"'^ at 
2 ^ z < 4.5. Since our analysis assumes a dust cor- 
rection consistent with these authors at high redshifts, 
but is based on more recent data, it thus seems that the 
SFR density peak suggested by their fit may be an arti- 
fact of the scarce high-redshift data in their compilation, 
which may be affected by cosmic variance and is not uni- 
formly complet e. For example, one of the z 6 points 
that driv es the IHopkins fc Beacorrj |2006) fit is the es- 
timate of iBunker et al.l ()2004[ ). which is only complete 
to O.li* and is based on an extremely small HUDF 11 
arcmin^ exposure. W e instead consider the analysis of 
iBouwens et al.l ()2007( ). which includes the HUDF data 
as a subset and yields a higher SFR density, and con- 
sistently integrate the LF down to zero luminosity. Al- 
ternatively, present surveys may be missing a significant 
UV luminosity density from very faint galaxies. 

It is immediately clear from the scatter in panels B and 
C of Figure [T] that the total UV luminosity density ex- 
trapolated from the measured LF should be interpreted 
with caution. In fact, the dispersion between different 
points at fixed redshift is generally larger than the calcu- 
lated error bars, indicating that these are unlikely to be 
uniformly reliable, a situation which is particularly man- 
ifest at z ^ 4. There are several reasons why this may be 
the case, including extrapolation to fainter magnitudes 
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than probed by individual surveys, cosmic variance aris- 
ing from large-scale structure, and parameters (perhaps 
inaccurately) held fixed in some fits. 

A number of previous studies of the LBG UV LF have 
also found little evid ence for a decHne o f the SFR density 
beyond z 3 (e.g., Steidel et al.l 119991 : iGiavalisco et al.l 
I200I lYoshida et al.ll2006D . This finding has in addition 
been corrobor ated by measures based on photometric 
redshifts (e.g.. [Thompson et al.|[200lt iThompsonI 120031) 
and is also in qualitative agreement with theoretical mod- 
els th at predict a SFR history peaking at higher redshift 
(e.g., ISpringel fc HernquistI l2003t iHernquist fc Springell 

mm . 

If only because the SFR density is expected to rise con- 
tinuously on physical gr ounds, it must event ually decline 
toward high redshifts. iBouwens et al.l (|2007f ) in fact find 
evidence for such a decline toward z = 6 on the basis 
of evolving dust obscuration suggested by observed /3- 
values at this redshift (e.g.. iStanwav et all [2005). We 
simply contend here that neither the present Lya for- 
est data or the recent UV LF compiled here, especially 
when considered together with their mutual scatter after 
extrapolation down to zero luminosity, show convincing 
evidence for the often-assumed peak in SFR density near 
z ~ 2 — 3. The requirement that the Universe be reion- 
ized by z = 6 also supports a SFR p eaking significantly 
earlier (|Faucher-Giguere et al.l[2008ah . 

5. COMPARISON WITH PREVIOUS WORK 

Similar conclusions have been reached in previous stud- 
ies of the UV background. , Bolt pn et al. (_2005) , in par- 
tic ular, inferred F fro m the Lya opacity measurement 
of ISchave et al.l (|2003l ) and also found its evolution to 
be consistent with being constant at 2 < z < 4. By 
comparing with the estimated quasar contribution, they 
also found evidence for a stellar-dominated UV back- 
ground at all redshifts. In their analysis using the 



iMiralda-Escude et al.l (|2000D PDF, iBecker et al] (|2007f ) 
also derived a flat F over this redshift range. Our re- 
sults extend beyond previous analyses in highlighting 
that common assumptions regarding the star formation 
history fall short of providing for the ionizing rate of the 
forest at z > 3. 

Measurements based on the proximity effect (e.g., 
IScott et al]|2000D have tended to yield F values higher by 
a factor of ~ 3. However, the overdense regions in which 
quasars reside are likely to bias these measurements high 
(|Faucher-Giguere et ai]|2008b[) . 

6. REIONIZATION 

The decline of the quasar LF and the increasing domi- 
nance of stellar emission to the high-redshift z > 3 ioniz- 
ing background make a compelling case that the Universe 
was reionized by stars. This gives credibility to analyt- 
ical and numerical calculations of hydrog en reionization 
that make this assumption (e.g., Furlanc tto et al.l l2004t 
iMcQuinn et~aI1 l2007bl : IZahn et al.ll2007f l. This is en- 
couraging news for upcoming observational probes of the 
epoch of reionization, such as redshift ed 21-cm emission 
and high-redshift Lya emit ters (e.g., IZaldarriaga et al.l 
l2004HMcQuinn et al.]l2007ar ). whose detailed interpreta- 
tion will rely on our understanding of the morphology of 
reionization and its origin. 
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